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1 Abstract 

A numerical study of a jet in crossflow with ground effect is conducted using OVERFLOW with 
dual time-steppig and low Mach number preconditioning. The results of the numerical study are 
compared to an experiment to show that the numerical methods are capable of capturing the dom 
inant features of the flow field as well as the unsteadiness associated with the ground vortex. 


2 Introduction 

The study of a jet impinging on a ground in crossflow has been the subject of both experimental 
and numerical 12 ’ 31 studies in the past due to interest in V/STOL vehicle behavior in near-hover 
conditions. Of interest are issues such as the extent, location, strength, and frequency of the 
ground vortex, which are crucial for the safety of the ground crew and for assessing hot gas inges- 
tion effects on aircraft engine operation. 

The flowfield associated with the present test case has several important features. When the jet 
impinges on the ground, the air flows radially outward from the point of impact. A ground vortex 
is created in the shape of a horseshoe due to the presence of a crossflow. Numerical studies of a jet 
in crossflow [2] and V/STOL aircraft [3> 4 ’ 5] have been successfully performed. In the numerical 
investigation of Ref. 2 , the flow did not exhibit unsteady behavior. In order to make it “go 
unsteady,” the jet needed to be pulsed. In the studies of Refs. 3, 4, and 5, experimental verifica- 
tions of the ground vortex predicitons were not available. In the current paper, the focus is on 
unsteady simulations of the jet in crossflow with ground effect for comparison with the experi- 
ments of Cimbala et. al. [11 
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The unsteady jet computations are obtained using the OVERFLOW code along with precondi- 
tioned dual-time scheme.™ 1 1 • Preliminary results shown in this abstract demonstrate that the 
gross features of the flowfield are well represented and the unsteady frequency associated w,th the 
Sound vortex agrees well with the experiment. The final version of the paper with tnclude gnd 
refinement studies as well as detailed parametric studies varying the crossflow velocity, jet veloc- 

ity, jet diameter, and jet height. 

The following section describes the numerical method used for the simulation of the jet. Sec. 
4describes the grid system and the problem configuration. In Sec.5, a numerical solution o J 
in Crossflow and ground effect is shown and the details of the solution are discussed. 

3 Technical Approach 

The OVERFLOW^ Navier-Stokes flow solver is used to perform the present computations. The 
diagonalized approximate factorization algorithm of Ref. 8 is used in conduction with dual time- 
stepping to get time-accurate results. A low Mach number preconditioner based on Refs. 9, and 10 
is used to capture the Mach 0.013 crossflow accurately. 

The implementation of a low Mach numbet precondttioner for unsteady flows in OVERFLOW is 
discussed in detail in Ref. 1 1 . For completeness, a summary of the formulation follows. 

3 1 Low Mach Number Preconditioning and Dual time-stepping 

The preconditioned dual-time stepping method in OVERFLOW is based on the dual toe-step- 
ping methods presented in references [7][10j. The Navier-Stokes equations can be written in con- 

servative form as 

de + dQ + 3£ + a£ w t(e) 

Pdx dt dx d y dz 

where q = L pu pv P w / is the vector of unknowns, E, F, and G are the inviscid fluxes and «fl) 
represents the viscous terms. The term r p is the preconditioning matrix. The artificial time term 

is introduced to the governing equations in order to provide a relaxation (sub-iteration) proce- 
3t 

dure between physical time steps. 

Termed dual time-stepping, the subiteration process is important for ensuring efficient conver- 
gence of the time dependent solutions at each physical time step. In combination with precon i- 
fioning it also allows the use of larger time steps to efficiently capture the low frequency 
“puffing” of the ground vortex. The preconditioning formulation is also important for insunng 
accuracy at the low Mach number conditions of the crossflow. 

In generalized coordinates, the Navier-Stokes equations are discretized with first order accurate 
Eufer implicit discretization for the artificial time term, second order backwards dlffe ^ ce ^ S . 
cretization for the physical time terms and central difference discretization for the spatial terms to 

obtain 
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, . , _ ,-i„ , _ r+ 3Ar r r = (a transformation matrix between conservative and 

wnere^ - j si ,& p - V 2Ar v 

primitive variables). A, B, and C are the flux Jacobians and 

= -A xT ^- 4 ^ 6 " - - 5 / + S / + 5 C 6 ‘ " Z - <e) ) ' 

The variable , is the time step counter while the variable k is the sub-iteration counter. Finally, 

40 - A diagonalized approximate factorization algorithm^] is used for the solution o 

equation When converged in artificial time, th.s method is formally second order accurate. 

The spalart-Allmaras turbulence model is used for the present computations. 

4 Problem Configuration 

The configuration is designed to match the conditions in the experiment of Ref. 1 . The experiment 
« conducted in a test sfctton where the top and a bottom wtdls « 21 M -be > 

inserted through the top wall and is 1 2 inches in length leaving a 9 inch gap before the jet hits the 
grmmd The jet exits a 3 inch diameter tube a. a speed of 150 ft/s or Mach 0 13. The crossflow 
^ed is 15 ft/s which is equivalent to a Mach number of 0.013. The Reynolds number based on 

the crossflow speed is 80000. 

The OVERFLOW code solves the Reynolds averaged Navier-Stokes equations on structured 
grids For the purpose of the numerical simulation, the jet profile was set to the idealized slug flow 
at the too of the grid A boundary layer was then allowed to develop in the jet tube and this was 
considered adequate for the pu^os/of obtaining a valid velocity profile at the exit of the jet tube. 

4.1 Grid system , 

Unlike a full aircraft simulation, the single jet has the advantage that there is no complex geome- 

hy'invoWed. Thus, a structured gnd in a single cylindrical block, shown m 
the Present numerical simulation. The structured gnd has the dimensions of 99(radial)x5 l(vert 
cal)x5 1 (circumferential) for a total of approximately 250,000 grid points The jet is placed at the 
center of the mesh This same jet and the wall of the jet tube are depicted in a cross section of the 
grid(99x51) shown in Fig. 2 to illustrate the problem and the boundanes. This cross section shows 
half of an XY cut through the grid such that the center of the cylindrical mesh corresponds to the 
right hand side of the figure. The incoming flow enters through the left boundary in ig. 
jet tube and the jet are on the right. 


5 Results . 

The typical flow field corresponding to a jet impinging on a ground is expected to consist of the jet 
hittinfthe ground and splaying in all directions. However, due to the crossflow P us ^ ^ 
ward moving part of the splayed jet backwards, a ground vortex develops in front of the je as 
shown in Fig P 3 . This ground vortex is blown backwards by the crossflow thus wrappmg Use 
around the j!t in a horseshoe shape when viewed from above. A particle trace of the jet and 
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crossflow is used to show this horseshoe shape in Fig. 4 . As indicated in Ref. 1 , the puffing of 
this ground vortex becomes the primary source of the unsteadiness. 



Figure 1: Oblique view of the grid 
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A common method to verify if a numerical solution has developed unsteadiness is to plot the 
residual of the discretized equations as a function of the time step. For steady flows, the residua 
tends to drop several orders of magnitude and then settle to a single value. In ideal cases, the 
residual is expected to go to machine zero. Unsteadiness in the flow field is clearly indicated y a 
residual that instead of settling to a single value, continues to increase and decrease. If the residu 
pattern has a specific frequency, it is most likely associated with some unsteady feature of the ful- 
filled. The residual history for the present test case is shown in Fig. 5 . The repeating pattern is 
seen in this history and corresponds to a frequency of approximately 1.4 Hz. The observed fre- 
quency in the experiments of Cimbala et. al. corresponded to between 1 and 8Hz with the fre- 
quency for the case presented here being approximately 3.5Hz. 


6 Summary 

The purpose of this paper is to numerically simulate this three dimensional unsteady flow in order 
to verify that the features in the flow field can be obtained numerically. Furthermore, computa- 
tions will be performed for several conditions in order to compare the unsteady frequency in the 
numerical solution to the frequency in the experiment of Cimbala et. al. 



Figure 3: A vector plot of the jet and the ground vortex 
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Figure 5: Residual history of the unsteady solution 
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